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Abstract 

Combining with our newest CCD times of light minimum of EM Cygni, all 45 available times of light 
minimum including 7 data with large scatters are compiled and the updated 0-C analysis is made. The best- 
fit for the 0-C diagram of EM Cygni is a quadratic-plus-sinusoidal fit. The secular orbital period decrease 
rate -2.5(±0.3) x 10~''s means that magnetic braking effect with a rate of mass loss via stellar wind, 
2.3 X 10~^^ MQyr~^ , is needed for explaining the observed orbital period decrease. Moreover, for explaining 
the significant cyclical period change with a period of ~ 17.74(±0.01)3'r shown in the 0-C diagram, magnetic 
activity cycles and light travel-time effect are discussed in detail. The 0-C diagram of EM Cygni cannot 
totally rule the possibility of multi-periodic modulation out due to the gaps presented after 25000 cycles. 
Based on the hypothesis of a K-type third star in literature, light trave-time effect may be a more plausible 
explanation. However, the low orbital inclination of the third body (~ 7°. 4) suggests that the hypothetic 
K-type third star may be captured by EM Cygni. But assuming the spectral contamination from a block 
of circumbinary material instead of a K-type third star, the third star may be a brown dwarf in case of the 
coplanar orbit with parent binary. 

Stars: cataclysmic variables; Stars : binaries : eclipsing; Stars : individual (EM Cygni) 



1 Introduction 



EM Cygni is a typical dwarf novae and classified as the Z Cam type subclass, since the usually regular erup- 
tion cycle is occasionally interrupted by irregular fluctuations of brightness with a low-amplitude ~ 0"'.4 
( Downes & Shara . 19931). Although extensive photometric and spectroscopic studies have shown that EM 
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Cygni is both an eclipsing binary and a double-lined spectroscopic binary (e.g. lMumford & Krzeminskil . ll969l : 



Kraft 



1964]) . its accurately absolute parameters are s till a problern due to t he contamination f rom a K2-5V third 



star a 



ong the line of sight, which is discovered by 



North et al. 



(119811) suggested its mass ratio is larger than unity. But 



(l2000h. 



North et al. 



Robinson 



( 1974 ) and 



Stover et al. 



(|2000l) pointed out that the previous works 



never considered the spectral contamination from a K-type third star, which leads to the radial velocity semi- 
amplitude for the mass donor star is increased from 135 kms'^ to 202 kms~^, and derived the mass ratio of EM 
Cygni is le ss than unity, which solved a problem that the mass transfer in EM Cygni can cause dynamically 



instability dStover et al. 



19811) . 



AH orbital period analysis of EM Cygni have indicated that its or bital period probably is decre asin g (IMumford 



980 



Beuermann & PakuU L 



1984 



Csizmadia et al. 



2008). But 



Beuermann & PakuU 



(Il984h and 



Csizmadia et al. 



(120081) suggested that it s decrease rate was small and not significant. In the updated 0-C diagram derived by 



Csizmadia et al. 



(|2008|) . a secular orbital period decre ase for EM Cygni c annot be significantly obtained due 



North et al. 



(|2000l) has discovered a K-type third star 



to the large scatters of the new 7 eclipse times. Since 
contributes approximately 16% of the light from the system, a new 0-C analysis for EM Cygni is important to 
verify the previous conclusion. 

Li this paper, 45 available times of light minimum from 1962 to 2008 are presented in Sect. 2, including a new 
data from our observation. Sect. 3 deals with the 0-C analysis for EM Cygni. Finally, the discussions of the 
possible mechanisms for orbital period change are made in Sect. 4. 



2 Observation of times of light minimum 

A new time of light minimum is obtained from our CCD photometric observations with the PI Vers Array 1 300B 
CCD camera attached to the 2.4-m RC telescope at Gao Meigu (GMG) observational base of Yunnan Observa- 
tory in China. It was carried out on November 25, 2008 in white light at 10-s exposure time. Two nearby stars 
which have the similar brightness in the same viewing field of telescope are chosen as the comparison star and 
the check star, respectively. All images were reduced by using PHOT (measure magnitudes for a list of stars) of 



the aperture photometry packa ge of IRA F. Since a small hump locat es in the center o: 
to the light curves obtained by 



Mumford & Krzeminski 



and 



Csizmadia et al. 



' eclips e, which is similar 
(120081), we removed this 



hump firstly and then used a parabolic fitting method to derive a new CCD time of light minimum. Based on 
the uncertainty of the parabolic fit, the error of our mid-eclipsing timing can be derived as ~ 0''.000072. 



The large uncertainties in the mid-eclipse timings for the CVs with low inclination suggest that the times 
of light minimum of EM Cygni in literature o btained by using photoelectr ic method without given errors, 
should have lower accuracy. Considering that 



Beuermann & PakuU 



(|1984() have derived the error of their 



photoelectric data as 0'^.0006, we arbitrarily adopted the equal errors of 0*^.0006 for the three photoelectric data 
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observed in the si milar observation conditio ns (IMumford I . ll980l) . and the larger errors of O'^.OOl for the previous 
data. In addition, 



Csizmadia et al. 



(|2008|) reported 8 new eclipse times including an amateur eclipse data in 
Heliocentric Julian Ephemeris Dates (HJED), which correspond to Ephemeris Time (ET), and all other data 
are in Heliocentric Julian Dates (HJD), which correspond to coordinated Universal time (UTC). However, for 
avoiding spurious trends in the orbital period analysis, we converted all 45 the eclipse timings in many different 
time systems to Barycentri c Julian Dynarn ical Date (BJDD), which correspond to Terrestrial Dynamical Time 
(TDT), by using t he cod e ( Stumpffi 1980 ). Since the recent 7 eclipse times without given errors observed by 



Csizmadia et al. 



(120081) present larger scatters (~ 0'^.008) than that of the previous data and the 0-C diagram in 
their paper never show the completed scatters of the 7 data points, we supposed an uniform error of 0''.003 for 
the 7 eclipse times. All 45 available times of light minimum covering about half a century are listed in Table 1 . 



3 Analysis of orbital period change 



For matching with the BJDD data we applied, the ephemeris in HJD derived by iBeuermann & PakuU I (|1984|) 
should be converted to BJDD, 

T„,in = BJDD 2437882.8607(3) + 0^.290909155(25)£, (1) 

which is used to calculate the 0-C values of EM Cygni. After linear revision, the new epochs and average 
orbital period were derived as, 

Tmin = BJDD 2437882.86141(18) + 0'^2909090099(33)£, (2) 

with variance cri = 2'^.8 x 10"^. Based on Eq 2, the calculated 0-C values are listed in column 6 of Table 1. 
Since the 0-C diagram shown in Fig. 1 presents a possible cyclical period change and all the previous orbital 
period analyses implied that a secular change in orbital period may be presented in EM Cygni, we attempted 
to use a quadric -plus-sinusoidal ephemeris to fit the data. The weights of data points are calculated from the 
errors of the mid-eclipse times. The least-square solution leads to 

(O - C)i = -l'^.90(+0^.34) X 10"^ + 2'^.64(+0'^.29) x lO '^E 
-3''.69(+0^.44) X IQ-^^E^ + l'^.83(+0^.28) x 10^^ 
sin[0°.01617(+0°.00003)£ + 353°.5(+8°.8)], (3) 

with variance cr2 = 2''.0 x 10"^, which is almost half of cti. The resulting 0-C diagram is sh own in Fig. 1 . The 



F-test used to assess the significance of quadratic and sinusoidal terms in Eq. 3 is suitable (jPringle L 119751) . as 
long as the form of parameter A is corrected to be 

A = — ^ (4) 

(r/(n-6)' 
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where n is the number of data. Thus, a calculation gives F(4, 39) = 8.4, which indicates that the fit is significant 
well above 99% level. However, since the small quadric term in Eq. 3 implies that the quadric ingredient in 
the 0-C diagram of EM Cygni may be not significant, a possible linear-plus-sinusoidal ephemeris for the 0-C 
diagram of EM Cygni is estimated to be 



(O - C)i = -8^.14(+2^.78) X 10"^ + 2^.07(+0''.44) x 10'^ E 

2^'.06(+0''.26) X 10"^ sin[0°.018634(+0°.000007)£ + 273°.4(+8°.6)], 



(5) 



with variance = 2'^.5 x 10 which is larger than ctj- In order to describe the significant level of the quadric 
term in Eq. 3, the parameter A of F-test describing by a similar formula as Eq. 4, 



A = 



(crj 



crl) 



CT\l{n-6r 



(6) 



indicates F(l,39) = 18.3. This means that the orbital period decrease in EM Cygni is significant well above 
99% level. In addition, the reduced value for Eq. 3 is calculated to be 1 .29 and the time span of the whole 
data set is more than two complete periods. Thus, this cyclical period variations should be significant. After 
removing the quadric element in the top panel of Fig. 1, the cyclical period changes are shown clearly in the 
middle panel of Fig. 1. Although the F- tests and the reduced all suggest that a significant cyclical period 
change can be detected in the 0-C diagram of EM Cygni, the gaps after 25000 cycles caused by the sparse data 
may mask a non-sinusoidal modulation. Therefore, the cyclical period change presented in the 0-C diagram 
of EM Cygni may be a quasi-period change with a period of ~ 17.74(+0.01)jr. According to Eq. 3, the orbital 
period decrease rate is calcul ated to be --2.5(+0.3) x 10~^^ j' s~^, which i s nearly two orders of magnitude larger 
than the previous works (e.g. iMumford L 1 198d : 



Csizmadia et al. 



2008) 



4 Discussion 



4.1 Secular orbital period decrease 

The significant orbital period decrease shown in Fig. 1 implies a secular variation of orbital angular moment 
in EM Cygni . The i ipdated spectroscop ic studies indicated that the mass ratio of EM Cygni is less than unity 



(lNorthetal.Ll2000l : 



Welsh et al. 



20071) . This means that the orbital period increase in EM Cygni is expected 



due to the mass transfer from the mass donor star to the white dwarf. Therefore, an available mechanism 
causing the observed orbital period decrease in EM Cygni should be considered. 

In order to explain the observed period decrease, the decrease rate of orbital period of EM Cygni can be 
described as 



P — P„,l -\- Pgr + Pmt + P, 



(V) 
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Figure 1: The {O - C)\ values of EM Cygni are fitted with the elements displayed in Eq. 3. The quadric and 
sinusoid curves are plotted in the dash and solid lines, respectively, in the top panel. After removing the quadric 
element, the {O - C)2 is plotted in the middle panel. The solid circles and diamonds d enote the data in li t eratur e 



and in our observation, respectively. The 7 eclipse times are plotted by open circles (|Csizmadia et al. 
The residuals and their linear fitted solid line are presented in the bottom panel. 



2008 ). 
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Figure 2: The mass and separation of the third-body in EM Cygni depending on its orbital inclination /' are 
plotted in the left and right panels, respectively. The solid star and diamond in both panels refer to a K-type 
star and a brown dwarf, respectively. 



where Pmu Pgr, Pmt and P„,i, are the orbital period change rates deduced by mass loss, gravitation radiation, mass 
transfer and magnetic braking, respectively . Since the orbital period of EM Cygni is long enough {P ~ 6.9/7), 



the gravitation radiation can be neglected (IVerbunt & Zwaan 
case (i.e. P,„/ = 0), then the Eq. 7 can be described as 



198l|) . If we only discussed the conservative 



P — Pmt "t" 



mb- 



(8) 



Although the low orbital inclination of EM Cygni implies that the bright spot at the outer edge of the disk caused 



by the mass transfer may be invisible in eclipsing light curves (IMumford & Krzeminski 1 . 1 19691 : 



1980 



Stover et al. 



Mumford 



19811) . EM Cygni classified as a Z C am subtype CV should have a high mass transfer 



rate, which is almost equal to the critical mass transfer rate dWarner I . ll995|) . By using the experiential formula 



dShafter & Wheeler . 



198 



Mcrit(Mi,P) 



2.8 X 10^^(1 +0.11^/^^1)3-^^ 

x[0.5 - 0.227/0^(0. ll/'/Mi)]'° V'-^3g5-\ 



(9) 



where P and M i are the orbital perio d and mass of white dwarf, respectively, and the updated system parameters 
of EM Cygni dWelsh et al. 1 . 120071) . the critical mass transfer rate Merit for EM Cygni is calculated to be ~ 
1.2 X \0~^Moyr~^. In addition, the accretion rate 5 x 10'^^ Moyr'^ from the disk to the white dwarf derived by 



Winter & Sion 



(|2003|) can be neglected because it is nearly three orders of magnitude lower than the estimated 
Merit of EM Cygni. Thus the mass transfer rate M^t ~ M^rit can produce the orbital period change with a rate of 



ML. 



8.6 X 10-^2^ s~\ 



A calculation withEq. 8 leads to P,nh = -3.3 x 10 s s . According to 



Tout & Hall 



(1199 1|) . the magnetic braking effect can be calculated by using 

Pmb 



Ra iMi + Ml ■ 

2{ — f— -M,nb, 

^ a ' M1M2 



(10) 



where R^, a and M,,,^ are the Alfven radius of secondary, the separation of the binary and the stellar wind mass 
loss rate, respec tively. Comb i ning t he orbital inclination i = 67° and the orbital elements asini = 1.52 x lO^'cm 



for EM Cygni ([North et al. L 120001) . we estimated the separation of binary is about 231Rq. Then R^M,nb is 



derived as 5.1 x 10 ^R^Moyr ^ Assumption that the Alfven radius of the K-type secondary is the same as that 
of the Sun(i.e. ~ 15i?o), the mass loss rate caused by the stellar wind is estimated to be 2.3 x lO'^^M^yr'^ . 



4.2 Quasi-period change 



Hall 



(|l989l ) has stated that the binaries with secondaries of later than F5 have cyclical period changes because 



of a solar-type magnetic activity cycle in the convective shell. As for EM Cygni, a K3-type secondary implies 



that a solar-type magnetic activity c ycle in its convect i ve she! 



0-C diagram should be considered (|Applegate . 



19921 : 



Lanza 



'or ex plaining the observed oscillations in the 



20061). Although a calculation for EM Cygni 



has suggested that Applegate's mechanism cannot interpret the observed cyclical period variation, the fact 
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that Applegate's model can be ruled out as a viable explanation does not mean that the hypothesis of period 
modulations driven by magnetic activity cycles in the secondary has to be discarded (ILanza & Rodono L 12002 : 



Lanza . 



2005Ll2006|) . Based on the gaps shown in Fig. 1, a magnetic activity cycles may be a well explanation 



for the quasi-period change in EM Cygni. However, 



Lanza 



(120061) also pointed out that the new improved 



model they proposed may be faced with a energy balance problem in the case of the secondary components of 
CVs due to their fast r otations. Mo reover, the updated plot of mass ratio vs. secondary component's spectral 
type as the diagram in 



Hall 



(|1989() for almost all t ype of binaries has d efinitely rules the magnetic activity 
out as unique cause for a cyclical period variation (ILiao & Oian 1 . 12010|) . Considering the discovery of the 
spectrum from a K-type third star (INorth et al. L |2000() . it is more plausible that the light travel-time effect 
caused by a perturb ations from a tertiary component is used to explain the observed cyclical variation. Since 



North et al. 



(I2OO0I) has detected the spectrum of a K-type third star in EM Cygni, the hypothesis of a strictly 
sinusoidal period change in its 0-C diagram is expected. The sinusoidal fit presented in Fig. 1 suggests that 
the orbital eccentricity of the third star should be zero, and the amplitude of sinusoidal curve determines the 
projected distance 0.32(+0.05)Ai7 from the third star to the mass center of the triple system. Then, according 
to the Third Kepler Law, the mass function of the third component /(ms, i) can be calculated to be 1.0(±0.5) x 
IO'^Mq. Upon that, it is worth y discussing v yhethe r the third body resulting in the cyclical period changes is 



the K2-5V third star derived by 



North et al. 



(l2000h or not. In the followin g calculations , a com bined mass of 



1.0(+0.09)Mo + 0.77(+0.08)Mo for the eclipsing pair of EM Cygni is used dWelsh et al. I . l2007h . 



4.2.1 a K-type third star 



Based on the calibration of MK spectral types, the mass of a K2-5V star can be estimated to be about 0.67Mo, 
which implies that the orbital inclination of the third body is very low {1 ~ 7°. 4). And then the distance 
from the mass center of system is derived as ~ 6.6AU. In this case, the triple system is odd and the third 
body may be originated from a straggler, whic h is captured by EM Cygni in a long time ago and now has the 
similar systemic velocity to the parent binary ([North et al. . 2000). Moreover, this geom etric structure of th e 



2000). 



triple system supports the prediction that the third star and the CV is at s i milar distances (iNorth et al. L 
Assuming the distances of EM Cygni is ~ 320pc (|Beuermann & PakuU L Il984h . the separation from the third 
star to the parent binary ~ 9.1AU projected on the sky can be estimated as ~ 0.03 , which is far beyond the 
resolution of the most advanced ground-based telescope. Therefore, this third star is hard to be separated from 
the parent binary. If this is the truth, then EM Cygni may be an important object for studying the dynamics of 
capturing a third body, and it is very interested to focus on the future evolution state of the K-type third star 
and the CV. 
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4.2.2 a brown dwarf as the third body 



Since a K-type third star in EM Cygni suggests an unusual triple system, we can discuss another possible 
case by assuming the spectral contamination for EM Cygni is not from a K-type third star but from a block of 
circumbinary material along the line of sight. Thus, the observed cyclical period changes should be caused by 
a new third body, which may be a brown dwarf with a distance ~ S.OAU from the mass center of system shown 
in the right panel of Fig. 3, as long as i' > 76° (with a low possibility < 15%). This implies that a brown dwarf 
as the third body is almost coplanar with EM Cygni and can survive the previous evolution of the common 
envelope evolution of the parent binary. In this case, EM Cygni hiding a brown dwarf is an imp ortant systern 
for fu rther study of the formation and evolution of substellar objects as the dwarf nova Z Cha (IDai & Oian . 



200% . 



5 Conclusion 



A new 0-C diagram of EM Cygni by adding a new eclipse time with high precision has been made and 
detailed orbital period analysis implies that the best-fit of the 0-C diagram shown in Fig. 1 is a sinusoidal 
variation with a period of ~ 17.74(+0.01)jr superimposed on the secular or bital period decrease with a rate of 
2.5(+0.3) X 10"^^ 5 s'^ . Using the updated absolute parameters of EM Cygni ([Welsh et al. I . l2007() . we proposed 
that magnetic braking effect may be a possible explanation for t he observed secular orbital period decrease. 
According to the experiential formula (|Shafter & Wheeler I . Il986|) . a critical mass transfer rate of Z Cam-type 
dwarf nova EM Cygni is estimated to be ~ 1.2 x 10"^Moyr"'. Moreover, we first detected a cyclical period 
change in the 0-C diagram of EM Cygni. Both plausible mechanisms that magnetic activity cycles and light 
travel-time effect are discussed. Although the gaps shown in Fig. 1 indicate that the for mer mechanism cannot 



North et al. 



20001) suggests that the latter 



be totally ruled out, the spectrum of a K-type third star detected by 
one also a possible mechanism. Based on the parameters of the K-type third star, we derived an odd structure 
of the triple system for EM Cygni with an extremely low orbital inclination of the K-type third star (/ ~ 7°. 4), 
which implies that the K-type third star may have a different origin (is captured by EM Cygni in a certain 
evolution stage) with parent binary. However, if we assumed that the spectral contamination is not from a 
K-type third star, but instead from a block of circumbinary material, then the observed cyclical period change 
suggests that the third body may be a brown dwarf whilst it is nearly coplanar with EM Cygni. We noted that 
the derived error of our mid-eclipsing timing seems to be underestimated for EM Cygni. This means that the 
conclusions we have made are just temporary and the more observations are expected. Anyway, EM Cygni 
is a controversial object and deserves to further observing in a longer base line including photometries and 
spectroscopies with high precision. 



8 



acknowledgements 



This work was partly Supported by Special Foundation of President of The Chinese Academy of Sciences and West Light 
Foundation of The Chinese Academy of Sciences, Yunnan Natural Science Foundation (2008CD157), and Yunnan Natural 
Science Foundation (No. 2005A0059M) and Chinese Natural Science (No. 10573032, No. 10573013 and No. 10433030). 
CCD photometric observations of EM Cygni were obtained with the 2.4-m telescope at Yunnan Observatory. We thank 
the referee very much for the helpful comments and suggestions that helped to improve this paper greatly. 

References 

Applegate, J. H., 1992, ApJ, 385, 621 

Beuermann, K., & PakuU, M. W., 1984, A&A, 136, 250 

Csizmadia, Sz., Nagy, Zs., & Borkovits, T, et al., 2008, AN, 329, 39 

Dai, Z. B., & Qian, S. B., 2009, ApJ, 703, 109 

Downes, R. A., & Shara, M. M., 1993, PASP, 105, 127 

Hall, D. S., 1989, SSRv, 50, 219 

Kraft, R. P, 1964, ApJ, 139, 457 

Lanza, A. P., & Rodono, M., 2002, A&A, 390, 167 

Lanza, A. F, 2005, MNRAS, 364, 238 

Lanza, A. P., 2006, MNRAS, 369, 1779 

Liao, W. R, & Qian, S. B., 2010, MNRAS, in press 

Mumford, G. S. & Krzeminski, W., 1969, ApJS, 18, 429 

Mumford, G. S., 1980, AJ, 85, 748 

North, R. C, Marsh, T. R., & Moran, C. K. J., 2000, MNRAS, 313, 383 

Pringle, J. E., 1975, MNRAS 170, 633 

Robinson, E. L., 1974, ApJ, 193, 191 

Shafter, A. W., & Wheeler, J. C, 1986, ApJ, 305, 261 

Stover, R. J., Robinson, E. L., & Nather, R. E., 1981, ApJ, 248, 696 

Stumpff, R, 1980, A&AS 41, 1 

9 



Tout, C. A., & Hall, D. S., 1991, MNRAS 253, 9 
Verbunt, F. & Zwaan, C, 1981, A&A 100, 7 

Warner B., 1995, Cambr. Astrophys. Ser. 28, Cataclysmic Variable Stars. Cambridge Univ. Press, Cambridge 
Welsh, W. R, Froning, C. S., & Marsh, T. R., et al., 2007, PASP, 362, 241 
Winter, L., & Sion, E. M., 2003, ApJ, 582, 352 



Table 1: The 45 times of light minimum for the Z Cam-type dwarf nova EM Cygni. 



BJDD 

24()()()()()+ 


type 


error^ 


Method 


E (cycle) 


(0 - Cf Ref. 


37882.860400 


pri 


.001 


pe 





-.001040 (1) 


37883.732200 


pri 


.001 


pe 


3 


-.002000 (1) 


37906.713100 


pri 


.001 


pe 


82 


-.002890 (1) 


37911.660400 


pri 


.001 


pe 


99 


-.001020 (1) 


37936.677900 


pri 


.001 


pe 


185 


-.001680 (1) 


37966.641400 


pri 


.001 


pe 


288 


-.001860 (1) 


37968.677900 


pri 


.001 


pe 


295 


-.001660 (1) 


37996.604900 


pri 


.001 


pe 


391 


-.001950 (1) 


38174.933600 


pri 


.001 


pe 


1004 


-.000470 (1) 


38345.698500 


pri 


.001 


pe 


1591 


.000826 (1) 


38348.605900 


pri 


.001 


pe 


1601 


-.000843 (1) 


38496.970200 


pri 


.001 


pe 


2111 


-.000129 (1) 


38561.552400 


pri 


.001 


pe 


2333 


.000273 (1) 


38562.424300 


pri 


.001 


pe 


2336 


-.000546 (1) 


38624.388600 


pri 


.001 


pe 


2549 


.000134 (1) 


38674.715700 


pri 


.001 


pe 


2722 


-.000020 (1) 


38675.588400 


pri 


.001 


pe 


2725 


-.000050 (1) 


38676.751400 


pri 


.001 


pe 


2729 


-.000689 (1) 


38878.934400 


pri 


.001 


pe 


3424 


.000561 (1) 


38883.879600 


pri 


.001 


pe 


3441 


.000314 (1) 


39052.604400 


pri 


.001 


pe 


4021 


-.002120 (1) 


39054.642900 


pri 


.001 


pe 


4028 


.000029 (1) 


39230.933600 


pri 


.001 


pe 


4634 


-.000122 (1) 


39293.770500 


pri 


.001 


pe 


4850 


.000438 (1) 


39767.662600 


pri 


.001 


pe 


6479 


.001630 (1) 


39769.697300 


pri 


.001 


pe 


6486 


.000016 (1) 


40006.788800 


pri 


.001 


pe 


7301 


.000674 (1) 
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Table 1: Continued. 



BJDD 
2400000+ 


type 


error^ 


Method 


E (cycle) 


(0 - cY 


Ref. 


40007.953600 


pri 


.001 


pe 


7305 


.001850 


(1) 


40008.826500 


pri 


.001 


pe 


7308 


.002050 


(1) 


41980.605600 


pri 


.001 


pe 


14086 


-.000162 


(2) 


41982.643700 


pri 


.001 


pe 


14093 


.001620 


(2) 


42515.879400 


pri 


.001 


pe 


15926 


.001130 


(2) 


43776.677400 


pri 


.0006 


pe 


20260 


-.000566 


(2) 


43778.714900 


pri 


.0006 


pe 


20267 


.000535 


(2) 


43780.751100 


pri 


.0006 


pe 


20274 


.000436 


(2) 


45257.409400 


pri 


.0006+ 


pe 


25350 


.004570 


(3) 


50692.461300 


pri 


.001 


pe 


44033 


.003480 


(4) 


53650.426600 


pri 


.003 


* 


54201 


.005960 


(4) 


53709.187300 


pri 


.003 


ccd 


54403 


.003090 


(4) 


53984.393700 


pri 


.003 


cod 


55349 


.009520 


(4) 


53989.334200 


pri 


.003 


ccd 


55366 


.004630 


(4) 


53990.500600 


pri 


.003 


ccd 


55370 


.007310 


(4) 


53991.367500 


pri 


.003 


ccd 


55373 


.001530 


(4) 


53993.408700 


pri 


.003 


ccd 


55380 


.006370 


(4) 


54796.020231 


pri 


.000072 


ccd 


58139 


-.000072 


(5) 



Note. * the given uncertainty in literature. * tlie unknown 
pe and cc d means the pliotoelectric and CC D methods used 



Csizmadia et 



Mumfo^^ Krzeminski 
al. |j2008l) ;(5)This paper. 



19691); (2) 



Mumford 



method used by 
in observations, 



amateur astronomer. 
respective ly. Refer- 



19801 : O^ lBeuermann & Pakulll 1198411 : (4) 
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